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We address the question of the Lorentz nature of the effective interquark interac- 
' tion in QCD which leads to the formation of the QCD string between colour charges. 

o ; 

■ In particular, we start from a manifestly vectorial fundamental interaction mediated 
^ , by gluons and demonstrate that, as soon as chiral symmetry is broken spontaneously, 

the effective interquark interaction acquires a selfconsistently generated scalar part 

■ which is eventually responsible for the formation of the QCD string. We demonstrate 
r~| . this explicitly for a heavy-light quarkonium, using the approach of the Schwinger- 

^: 

' , Dyson-type equation and the quantum-mechanical Hamiltonian method of the QCD 

D 



string with quarks at the ends. 
PACS numbers: 12.38.Aw, 12.39.Ki, 12.39.Pn 



I. INTRODUCTION 

In this paper we discuss one of the long-standing problems in QCD — namely, the 
problem of the Lorentz nature of the long-range confining interquark interaction. QCD 
is believed to be a stringlike theory at large distances, that is the long-range interquark 
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■ interaction is expected to be generated by an extended object — the QCD string. Such a 
phenomenological picture appears rather successful in various studies of hadronic properties. 
An important step was made in the framework of the Vacuum Correlators Method (VCM) 
|l| in which the Lagrangian of the QCD string with quarks at the ends can be derived natu- 
rally starting from the fundamental QCD Lagrangian 2|. In addition, the nonperturbative 
spin-dependent forces in heavy and light quarkonia were found, following the formalism 
established in Ref. Thus, for the spin-orbit interaction, a celebrated representation, 

V^^^(^) = I M _ ^\ (l^ + l^) + 1 (ff_i _ f\ 1^ (1) 
4m? j \r dr r dr J 2mgmq V ^ "^Z r 9r ' 



was introduced in Ref. j^, where e{r) is the static confining potential, and a general relation 
(Gromes relation j2]) is valid, 

e' + VI - V2 = 0. (2) 



For a purely scalar interaction, one obtains at large r's: 



v; = -e', ^' = 0, 



(3) 



and this was demonstrated explicitly for the Gaussian approximation for the field correlators 
in Ref. j^. For the case of a vector confinement, for example, for the Coulomb potential, 
one would find: 

v; = o, v^>o, (4) 

and the coefficient at the spin-orbit term would have the opposite sign. Phenomenology 
of the heavy quarkonia spectrum favours the first possibility, Eq. ([3]), so that one has an 



:or heavy quarks. On the 
6|] for recent results and 



evidence that Nature prefers scalar interquark interaction, at least 
lattice, numerous data also support the first possibility (see Ref. 
the vast bibliography). In the meantime, any quantum-mechanical approach meets severe 
problems with the description of another celebrated phenomenon — spontaneous breaking 
of chiral symmetry, which is known to take place in QCD. A full quantum field theory 
based treatment has to be exploited for this purpose. An example of such a treatment, 
also based on the VCM, is given by the Schwinger-Dyson-type approach to heavy-light 



quarkonia suggested in Ref. [7], and the 



jorentz nature of confinement for heavy quarks was 



studied in this formalism in Refs. 7, [s], ol. On the other hand, for light quarks, vectorlike 
confining interaction would have resulted in the well-known Klein paradox and, thence, in 
problems with building the spectrum of hadrons. No evidence for such problems exists so 
far. In this paper we prove that, indeed, even for light quarks, if chiral symmetry is broken 
spontaneously, the effective interquark interaction acquires a selfconsistently generated scalar 
part. This result is quite general and holds regardless of the explicit form of the interquark 
kernel, suffices it is confining and thus leads to spontaneous breakin g of chiral symmetry. A 
link is established between the VCM and potential quark models 10|, lul] which we refer as to 
Generalised Nambu-Jona-Lasinio (GNJL) models, the latter being widely used for studies 
of low-energy phenomena in QCD. This is an important outcome of our work since the vast 
results obtained in the literature in the framework of such quark models are valid for our 
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situation as well — what we do is approaching the same problem from another side. Thus we 
demonstrate that, starting from the fundamental QCD Lagrangian and using the Gaussian 
approximation for the background field correlators, one can derive a Schwinger-Dyson-type 
equation for the heavy-light quarkonium which, at large interquark distances, reduces to 
a Diraclike equation with an effective interquark interaction which contains a dynamically 
generated scalar part, as a consequence of chiral symmetry breaking. At the same time, a 
Schrodingerlike equation with the Hamiltonian of the QCD string with quarks at the ends 
(in the form of the Salpeter equation) arises naturally from the same Schwinger-Dyson-type 
equation, if the scalar interaction dominates We conclude, therefore, that this is the 
dynamical scalar interaction responsible for the QCD string formation. 



II. SCHWINGER DYSON-TYPE EQUATION FOR A HEAVY LIGHT 

QUARKONIUM 

In this section we consider a heavy-light quarkonium consisting of a static antiquark and 
a quark, whose mass is unconstrained and can take any value (we shall be mostly interested 
in the case of the massless quark). Our starting point the heavy-light Greens function SqQ 
written in Euclidean space as j?] 

S,Q{x,y) = D^pDi^^DA^exp |-i J d'xF^l - j d^x^j^-td -%m- i)^| (5) 



^^\x)SQ{x,y\A)^{y), 

where S'q(x, ?/|A) is the propagator of the static antiquark placed at t 
analysis it is convenient to fix the modified Fock-Schwinger gauge [13 



le origin. For further 



xA{x^,x) = 0, Ai{xi,6) = 0, (6) 

which ensures that the gluonic field vanishes at the trajectory of the static particle. As a 
result, the static antiquark decouples from the system, it is Green's function being simply 

Sq{x, y\A) = Sq{x, y) = i^-^d{x, - y,)^'^'^^--^'^^ + ^^^^(?/4 - x^)e-^'^y-~^'^\ (7) 

It is easy now to perform integration of the the gluonic field in Eq. ^ to arrive at 

S^Q{x,y) = ^ / D^D^^exp |- / d^xLeff(^,^^)|^^(a;)5Q(a;,t/)^(y), (8) 



c 
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with Leff(V',V'^) being the effective Lagrangian of the hght quark moving in the field of the 
static antiquark source: 



2 „ 

where all a's and /3's are fundamental colour indices, and the irreducible correlators 
(y4^^^^^(xi) . . . of all orders enter. The first correlator, {A^'^), obviously van- 

ishes due to the gauge and Lorentz invariances of the vacuum. In what follows we assume 
the Gaussian dominance to take place in the QCD vacuum and thus we keep only the bilocal 
correlator 

= 2(A,)^(A,)Jir,.(x, y), (10) 

and neglect contributions of all higher correlators^. Then, using the relation (Aa)^(Aa)]' = 
— 2j^'^i3'^J considering the large-iV^ limit, we rewrite the effective light-quark 
Lagrangian in the form: 

L,^{^^^^)=^l{x){-id-im)r{x) + ]^ j d'y i^i{x)^,i^^{x)i^l{y)^,r{y)K,.{x,y), (11) 

which leads to the Schwinger-Dyson-type equation Q, 

(_.4_„„)5(,,,)_.J,.,M(.,.)S(..y) = ^H)(,_,), 

-iM{x,z) = K^i,{x,z)-ff,S{x,z)-f^, 



(12) 



for the colour trace of the light-quark Green's function S{x,y) = ■^{ip'^{x)4'p{y)) . 

In order to evaluate the quark kernel K^y{x,y) we notice that a celebrated property of 
the radial gauges (gauge ([6]) obviously belonging to this class) is a possibility to express the 



^This approximation leads to an exact Casimir scaling, that is the ratio of any two potentials between 
static sources in different representations of the colour group is given by the ratio of the Casimir operators 
evaluated for the given two representations [14]. The Casimir scaling was tested on the lattice [15] and 
it was found to manifest itself with a very high accuracy, which evidences a suppression of higher gluonic 
correlators as compared to the Gaussian one and thus justifies the approximation made above. 
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gluonic field A in terms of tlie field strength tensor F. For the gauge ([6]) such a relation 
reads: 

A1(x4,x) = / daXiF",(x4, ax) 



\ (13) 

A^{x4,x) = / axkFI^^{x4, ax)da, i = l,2,3, 
Jo 

so that the kernel K^j^^, can be expressed in terms of field strength correlator {F^^{x)F^p{y)), 
for which we use the parametrisation 

(C(^)^t(l/)) = j;^D{x - - 5,/.,) + A(i), (14) 

where the second term A^^^ is a full derivative and it does not contribute to confinement 
and therefore will not be considered below. The profile function D{x — y) decreases in all 
directions of the Euclidean space, and this decrease is governed by the gluonic correlation 
length Tg. Lattice simulations give rather small values of Tg ^ 0.2 0.3 fm [g], [igI, so that 
the profile D{x — y) has the support at close point y ^ x. The term proportional to D{x — y) 
in contributes to the area law with the string tension |l| 

/■oo 

dT / d\D{T,\). (15) 

Then, in view of Eqs. (fT3l) and (iMl) . the quark kernel K^y{x,y) = K^y{xi — yi,x,y) is 
found to be {t = x^ — y^) 

K44{T,x,y) = (xy) da d(3D{T,\ax - (3y\), 
Jo Jo 

Kii{T, X, y) = K4i{T, X, y) = 0, (16) 
Kik{T,x,y) = {{xy)6ik - yiXk) ada (3d(3D{T,\ax - (3y\). 



The Schwinger-Dyson-type equation f[T2|) is an essentially nonlinear equation. Its lin- 
earised form, with the Green's function S{x, z) substituted by the free-quark Green's func- 
tion 5'o(x, in the mass operator M{x, z), can be used if the quark is heavy. This was done 
in Refs. |8|, l9[], where the effective potential for the heavy quark interaction with the static 
antiquark source was built, including the nonperturbative spin-orbit interaction. The lead- 
i„g correction to this potenfa. due to the proper string dynamics was identified .n Ref. Q. 
It was noticed in Ref. [9|, however, that the given linearisation procedure is selfconsistent 
only if the product of the quark mass and the gluonic correlation length is large, mTg ^ 1. 
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In case mTg -C 1 the series of corrections to the leadin 



regime blows up and no conclusion 



concerning the dynamics of the system can be made |9|]. This procedure is useless therefore 
for the purposes of the present paper which is aimed to consideration of the light (massless) 
quark with its effective mass generation due to the phenomenon of spontaneous breaking of 
chiral symmetry. Thus we study this case using a different approach. 

Below we use two simplifications: i) we neglect the spatial part of the kernel Kij^ and ii) 
neglect corrections due to the finitness of the correlation length Tg. The first approximation 
utilises the fact that, although the spatial part of the kernel is important for the correct 
account of the QCD string rotation, it is not decisive for the Lorentz nature of confinement, 
yielding only unnecessary complications. The second approximation is justified in view of 
the results of the lattice simulations which give, as was mentioned before, quite small values 
of Tg. The latter simplification allows us to approximate the kernel ( IT6|) by an instantaneous 
kernel and thus, for the Fourier transform of K in time, to neglect its dependence on the 
energy, 

/•l /•! POO 

i^44(cu, X, = i^(cu, X, = = (xy) / da d(3 I dTD{T,\ax — (3y\). (17) 

JQ Jo J-oo 

Then the mass operator can be written as 

M(x, y) = (5(x4 - y^)M{x, y), M{x, y) = ^K{x, yhMx, y), (18) 

where, following Ref. [3], we introduced the quantity 

°° f duj 

H^,y)= ^nix)sign{n)ijl{y) = 2i / —S{uj,x,y)ji = 2iS{x4-y4,x,y)-f4\^^^y^, (19) 

n=—oo 

which is convenient for studies of the Lorentz nature of the interquark interaction. In 
addition, in the limit Tg — » 0, the profile function D{t, A) takes a singular form D{t, A) = 
2cr(5(r)5(A) (see also Ref. [18|| for the discussion of the singular limit of some stochastic 
model). For such a profile one finds readily for the kernel f|T7|l : 

K{x,y) = 2a{xy) [ da! d(3 6{\ax - /3y\). (20) 
Jo Jo 

Evaluation of the integrals in Eq. (I2U]) is trivial and gives 

K{x,y} = 2(jmin(|f|, \y\) = a{\x\ + \y\-\x-y\), (21) 
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if the vectors x and y are collinear (the kernel vanishes otherwise, as required by the delta- 
function in Eq. fl20|) ). The requirement of coUinearity of the vectors x and y ensures that 
the interaction between the light quark and the static antiquark is due to an infinitely thin 
string — a two-dimensional object embedded into the four-dimensional space. In order to 
proceed we simplify the form of the kernel (I2T1) and relax the constraint of collinearity. Thus 
we approximate the kernel as 

K{x,y) = a{\x\ + \y\ - \x-y\). (22) 
The quark kernel (!22l) possesses a number of important properties: 

• it allows us to pass over back, to Minkowski space — it will be used from now onward 
in this paper; 

• it admits a clear interpretation. Indeed, the kernel can be split into two parts: the 
local part —a\x — y\ which is responsible for the selfinteraction of the light quark, and 
the nonlocal part cr{\x\ + \y\) which describes the interaction of the light quark with 
the static source. Such a form of the kernel is a consequence of the gauge condition 
([6]) which decouples the static particle from the system and brings all the information 
about the antiquark to the kernel K{x, y); 

• it admits a natural generalisation from the linearly rising potential to the potential of a 
generic form V{r). In order to emphasise this important property we keep the potential 
as V{r) in all formulae below. Nevertheless, every time we need to specify the form of 
the potential, the linear confinement is understood, as the most phenomenologically 
justified candidate; 

• with the kernel (|22ll we establish a link between the VCM and the GNJL models for 
QCD with instantaneous quark kernels which have a long history in the literature 
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11| and which can be viewed as nonlocal divergence-free generalisations of the 



Nambu-Jone-Lasinio model [19(] (recent exhaustive studies of the mesonic spectrum 
in this model can be found in Ref. Below we employ the chiral angle approach 

which is widely used in such potential quark models. 

Combining Eqs. f|T2l) . (fTSl) . f|T9l) . and fl22|) together we arrive at the Schwinger-Dyson-type 
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d d 

-ho^ + - m ] S{t, x,y)- I £zM{x, z)S{t, z, y) = 5{t)5^^\x - y), (23) 



equation for the heavy-light quarkonium in the form: 

d _ a 

' dt dx 
where 

/duj 
^5(^,f,y)7o, (24) 

and the quark kernel is 

K{x, y) = V{\x\) + V{\y\) -Vi\x-y\). (25) 
In the next section we study the properties of this equation. 

III. THE LORENTZ NATURE OF CONFINEMENT 

Let us investigate the properties of the quantity A{x,y) and demonstrate the way it 
acquires the contribution with the matrix structure oc 70, since exactly this phenomenon 
constitutes spontaneous breaking of chiral symmetry. 

It was argued in Ref. [?] that Eq. (fT2!) admits linearisation via the substitution 

A{x,y)^-fo6^'\x-y) + ..., (26) 

where the ellipsis denotes subleading at large distances terms. Such a substitution was 
justified then by an explicit computation of this quantity using the spectrum of the resulting 
linearised Eq. (123|1 . In Appendix [A] we give some details of the derivation of Eq. (!26|) . taken 
from Ref. [3] • The scalar Lorentz nature of the effective interaction follows immediately from 
the form of Eq. (1231) and the matrix structure of A{x,y) given in Eq. (|26l) . We conclude 
therefore that Eq. ( l23l) does admit a solution given by the scalar interaction generated in a 
selfconsistent manner. Let us have an insight into the details of this selfconsistent generation 
of the scalar effective interquark interaction. 

The separation of the kernel into the local and nonlocal parts mentioned above hints us 
the way to proceed. Indeed, let us consider the local part of the kernel first and omit the 
nonlocal part. Then Eq. (|T2l) reduces to the Dyson equation for the light-quark propagator 

(70P0 - 7P- - S(p))5'(po,p) = 1, (27) 
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where the mass operator S(j)) does not depend on the energy due to the instantaneous 
nature of the interaction and can be evaluated as 

m = -^J ^4^iP- khoSih, kho. (28) 

Equations (1271) an d (125]) together lead to the selfconsistent nonhnear equation for the 
quark mass operator [lOj, 

^'^^--'/|^^'^--^>° -,„,„-,g!„-,(,V °- '''' 

If we parametrise now the mass operator in the form: 

S(p-) = [A, - m] + m[B, - p], (30) 

with Ap and Bp being two auxiliary functions, then Eq. (!29|) gives the selfconsistency con- 
ditions for such a parametrisation, 

Ap = m+^ J -^^V{p-k)smipk, Bp = p + ^ J {pk)V{p- k) cos (fk, (31) 

where the angle ipp — known as the chiral angle — is introduced to obey the condition 

Ap cos ipp = Bp sin ipp, (32) 

which, given the relations (13T!) . plays the role of the mass-gap equation for the chiral angle. 
Historically the chiral angle ipp is defined such that (pp{p = 0) = | and ipp{p oo) = 0. In 
Fig. [U we plot the chiral angle — solution to the mass-gap Eq. (!32l) for the linear confinement. 



The interested reader can find the details of the chiral angle formalism in Refs. 



3, 
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2l| 



A comprehensive analysis of the properties of 
various powerlike potentials is given in Ref. 



;he mass-gap equation and its solutions for 



22 



It is an easy task now to evaluate the quantity A{p,q), which is the double Fourier 
transform of A{x,y): 

m q) = 2tj ^S{uj,p, g-)7o = {27Tf6^'\p- q)Up, (33) 

where 

Up = p sin Lpp - {ap) cos v^p, /5 = 7o, a = 7o7- (34) 

Equation (IMI) gives the answer to the question on the Lorentz nature of confinement 
in the heavy-light quarkonium. Indeed, for low-lying states with the relative momentum 
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p being small, the cliiral angle (pp is close to tt/2, so that the matrix Up = (3 and this 
immediately leads one to Eq. (1261) . Notice that the contribution to A{x,y) proportional to 
the matrix 70 appeared entirely due to chiral symmetry breaking described in terms of the 
nontrivial chiral angle ifp — see Fig. [1] {(fp = for the massless quark and without chiral 
symmetry breaking^). This is the regime found in Ref. [3] and mentioned in the beginning of 
this section. In the next section we demonstrate that one needs exactly this regime to realise 
in order to be able to describe the quarkonium using the Salpeter equation. The opposite 
situation of the vanishing chiral angle, which realises for highly excited bound states, is 
discussed in detail in Ref. [23|. 



Obviously, inclusion of the spatial part of the kernel f|T6|) as well as relaxing other simpli- 
fying assumptions made in course of this section do not change the main conclusion of this 
section — namely that, as soon as chiral symmetry is broken spontaneously, an effective 
scalar interaction appears in a selfconsistent manner. 



^The situation is trivial for the heavy quark when chiral symmetry is broken explicitly. Indeed, in this 
case the chiral angle acquires the contribution arctan(TO/p) and thus ipp fa tt/2 for p <ti m. 
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IV. SCALAR CONFINEMENT AND THE SALPETER BOUND STATE 

EQUATION 



In the previous section we considered the approach to the heavy-hght quark-antiquark 
system based on the Schwinger-Dyson-type Eq. ( |T2l) . An alternative approach to (heavy- 
hght) quarkonia, based on the spinless Salpeter equation, 



[\/p^~-i-m^ + ar\il) = Eip, 



(35) 



is also celebrated in the literature (see, for example, Ref. 2J] and Appendix [B for the 
derivation of Eq. fl35l) in the^ormalism of the QCD string with quarks at the ends {2I, which 
is also based on the VCM p^]). The purpose of the present section is to demonstrate that 
Eq. fl35l) and its more sophisticated versions, like the Hamiltonian of the QCD string with 
quarks at the ends (see Appendix [B]), are consistent with the bound-state Eq. (HOj) under 
the assumption of the scalar confinement dominance in the effective interquark interaction. 

We return now to the full Eq. (l23l) and rewrite it in the form of the bound-state equation 
for the bispinor wave function \E', 



{dp + /3m)*(f) + (3 d^zM{x, z)^(i') = E^(f) 



(36) 

with the nonlocal part of the kernel included together with the local one. Notice that the 
heavy-light bound-state equation in the form Eq. ( l36i) can be derived independently in the 
formalism of the GNJL quark models {2^. 

Passing over to the momentum space and using the mass-gap equation in the form: 



EpUp = ap + l3m + 



(2^ 



V{p-k)Uk, 



one can rewrite Eq. (15^ as 



EpUp<^>{p) + \j 0^,V{p~ k)[Up + U,]^{k) = E^>{p), 



(37) 



(38) 



where the quantity Bp = Ap sin ipp + Bp cos ipp is the full quark dispersive law which substi- 
tutes the free dispersion \/fP + and which appears as a result of the quark selfinteraction. 
Equation (138|1 is subject to a Foldy-Wouthuysen transformation, which was built in a 



closed form in Ref. 



The corresponding operator is: 



T 



exp 



m{p) = Tp 



i){p) 



(39) 
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and the resulting Shrodingerlike equation, which stems from Eq. (138|) after the Foldy- 
Wouthuysen transformation with the operator fl39l) . reads: 

j^V{p- k) [CpCk + m{ak)SpSk\ m = Eij{p), (40) 

where we used the shorthand notations Cp = cos |(f — '^p) and 5*^ = sin |(f — v^p); are 
Pauh matrices, and p and A; are unity vectors for p and k, respectively. 

Let us consider Eq. fj40|) in the regime ipp ^ ^. Then Cp = 1, 5*^ = and the interaction 
term in Eq. fHU]) reduces to the plain potential ar, in coordinate space. Then the resulting 
equation reads: 

[Ep + ar]tp = EtP, (41) 

where Ep plays the role of the kinetic energy operator for the quark. This equation has the 
form of the Salpeter Eq. ( l35l) . Notice, however, that it is not always sufficient to keep the 
kinetic term for light quarks in the form of the free-particle energy \/fP + m^, as in Eq. fl35l) . 
Strictly speaking, the dispersive law of the light quark Ep is generated dynamically and has 
to be treated with care. This is especially important for small interquark momenta, where 
Ep can even become negative. As a result, the lowest states in the spectrum — the pions and 
the kaons — cannot be described using the simple Eq. and the like. These states are to 
be considered using either the full Bethe-Salpeter equation with the two-component mesonic 



wave function 
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2l| or in the framework of the full Schwinger-Dyson-type equation, similar 



to the heavy-light Eq. f|T2l) |25l . |26| . A progress in adapting the Salpeter equation based 
approach to description of lightest mesons was achieved in Ref. [27| in the framework of a 
matrix Hamiltonian technique. Apart from the aforementioned problem with the pions and 
kaons, the quark dispersive law Ep in Eq. (141 p can be substituted, with a good accuracy, by 
the free-quark energy, so that the Salpeter equation ( l35l) is readily reproduced. 

In addition, as seen from Eq. fHOj) . the simple Salpeter Eq. fl35|) has to meet certain 
problems for highly excited states as well. Indeed, the relative interquark momentum is large 
in excited mesons, so that the chiral angle vanishes asymptotically (see Fig. [1]). As a result, 
the effective interaction in Eq. (HUil becomes vectorial (see Eq. (HUI) with Cp = Sp = 1/ -\/2) 



2m- 



and it does not reduce to a plain potential anymore 

We conclude this chapter stating that, contrary to naive expectations (potential is added 
to the energy), the form of the Salpeter Eq. (!35l) does not suggest that the effective interquark 
interaction in the meson is vectorial. Furthermore, we demonstrate that this equation arises 
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naturally from the full Schwinger-Dyson-type equation for the heavy-light quarkonium un- 
der the assumption that chiral symmetry is broken (explicitly or spontaneously) and, as a 
result, the effective scalar interquark interaction is generated in a selfconsistent manner. 



V. CONCLUSIONS 



In this paper we address the problem of the Lorentz nature of confinement in QCD. We 
consider a heavy-light quark-antiquark system as a testground and exploit the Schwinger- 
Dyson-type equation derived for the Green's function of such a system using the VCM. 
We demonstrate explicitly that the stringlike picture of the interquark interaction at large 
distances (in the form of the Salpeter equation (135|1 ) appears due to chiral symmetry break- 
ing. In particular, we prove that the Salpeter equation fl35|) appears selfconsistently in the 
Schwinger-Dyson approach to the heavy-light quarkonium if chiral symmetry is broken, 
explicitly or spontaneously, and the effectively generated scalar potential dominates in the 
effective interquark interaction. This implies that the genuine Lorentz nature of the confin- 
ing interaction in this Salpeter equation (as well as in the Hamiltonian of the QCD string 
with quarks at the ends) is scalar. This is the main result of this work. This solves the prob- 
lem of the Klein paradox which is known to operate for systems with vectorial interaction 
growing with the distance. We conclude that there is no room for such a problem in QCD. 

The reported result is robust since it is only based on quite a general consideration and 
is stable across the whole variety of quark kernels. Furthermore, our conclusions acquire 
additional support from the fact that exactly the same bound-state equation for the heavy- 
light quarkonium can be derived independently in the framework of the GNJL quark models 
which have a long history in the literature and are known to give deep insight into physics 
of chiral symmetry breaking. 

Finally, let us mention that for light quarks and without spontaneous breaking of chiral 
symmetry one would have a vanishing chiral angle and, consequently, no effective scalar 
interaction. This situation is believed to realise in QCD above the temperature of the chiral 
symmetry restoration transition or for highly excited hadrons (see Ref. 28|] for a review). 
Properties of the interquark interaction in these situations deserves a special investigation 
and will be subject for future publications (see, for example, Ref. |29l|). 
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APPENDIX A: DERIVATION OF EQ. ([26]) 

Following the approach suggested in Ref. we assume (and justify this assumption 
a posteriori) that the Schwinger-Dyson Eq. fl23|) possesses a solution which gives for the 
quark mass operator (see Eq. (12^ ) a form described by a local scalar {U{x)) and local 
vector iV{x)) potential, 

M{x,y) = [U{x)+^oV{x)]S'^'\x-y). (Al) 
Then one can rewrite Eq. fl23!) in the form of a Dirac equation for the wave function ip{x), 

{ap + f3[m + U{x)] + V{x))iIj{x) = E^{x) (A2) 
or, in components {ip = Idp^Q^'^))^ 



^ + -G„-(E„ + m + [/-r)F„ = 
ar r 

^--F^ + {E^-m-U- V)Gn = 0. 
ar r 

One can use a simple trick to guess the matrix structure of the function A[x,y), 



(A3) 



M^^y)= ^ ^n{x)sign{n)iljl{y), (A4) 

n=—OD 

built with the help of the solutions to Eq. (]A2p . Indeed, according to its definition, A{x,y) 
can be naturally split into two parts, 

A(^)(f,f) = A^^\x,y) - A^Y\x,y), (A5) 

where ± stand for the summation over positive and negative eigenvalues, respectively. Also, 
for future convenience, we used the script {V). A similar decomposition is valid for the 
reversed sign of the vector interaction, 

A(-^) (f , y) = A^-^) (f , y) - A^-""^ (f , y) . ( A6) 
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Now we notice the following symmetry inherent to the system (1A3P : {V, En, k,,G, F) 
{—V, —En, —K, F, G) and find 




Hence one can rewrite Eq. flASP as 

A(V) = [aJ^) _ a(_-^)] + [aL-""^ - aL""^] = 70 J] (GG* - FF*) + 5k, (A8) 

where the correction 5A vanishes for = 0. Therefore, for a purely scalar confinement, 
the matrix structure of A{x,y) is, indeed, given by the matrix 70. In order to establish 
its spatial structure one can use the WKB calculation performed in Ref. j^]- We omit the 
lengthy calculation which can be found in Ref. p] and quote here the final result: 

A{x,y) ^7o n^, — Ko{ay^\x -y\)6{l -cosO^y), (A9) 

where Kq is the MacDonald function. It is easy to check that, 

AA(x,y) =70, (AlO) 



and, therefore, for \y\ ^ a\i-y\ ' ^i^^y) '^^^ t>e approximated by the three-dimensional 
delta-function peaked aX y = x. Thus we arrive at Eq. fl26l) . Moreover, for V ^ the same 
WKB method reproduces Eq. (lA9p and gives the decrease of the term 5A in Eq. (lASP at 
large distances, so that Eq. ( !26l) holds 

APPENDIX B: ROTATING QCD STRING AND THE SPINLESS SALPETER 

EQUATION 

In this appendix we give a brief derivation of Eq. ( l35l) in the formalism of the QCD 
string with quarks at the ends which is also derived in the framework of VCM. Following 
the method of Ref. {2], we start from the in- and out-states of the quark-antiquark meson, 

v|/(J'°"*)(x,i/|A) = ^,-(x)$(x,i/)vl/,(y), ^x,y) = Pexp(^tg dz^Ap'^y (Bl) 

and build its Green's function, 

G,,= {^f^'\x,y\A)^'^ff{x,y\Ay),,A = {TTS,{x,x\AMx,y)S,{y,y\AMy,x))A, (B2) 
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where Sg and Sq are the propagators of the quark and the antiquark, respectively, in the 
background gluonic field. Averaging over the background field can be done using the minimal 
area law assumption for the isolated Wilson loop, 



TrP exp ( ig (b dz^A 



C / / A 



exp (— o"5'n 



(B3) 



which is usually assumed for the stochastic QCD vacuum (see, for example, Ref. and is 
found on the lattice. Here S'min is the area of the minimal surface swept by the quark and 
antiquark trajectories. 



S'n.in = i dt [ dl3^{wW>Y - W 



2yjl2 



(B4) 



'0 ^0 

where, for the profile function of the string w^(t,P), we adopt the straight-line ansatz: 



W^{t,P)=Pxi^it) + il-P)x2, 



(B5) 



xi,2it) being the four-coordinates of the quarks at the ends of the string. We choose to 
consider the system in the laboratory frame and also to synchronise the quark times. 



Xio — X20 — t- 
The resulting Lagrangian of the string reads: 

i^str = -ar f d(3Jl-[n^{(3'xi + {l-(3)'x2)]\ 



(B6) 



(B7) 



where r = Xi — X2, n = r/r. This interaction Lagrangian is to be supplied by the quark 
kinetic terms —mi\ l 



x^ — m2 Y 1 — X2- Then, with the help of the auxiliary (einbein) 
field technique, used to get rid of the square roots in the kinetic (the einbeins fii^2 2^) and 
in the string term (the continuous einbein I'iP) I2I) in the Lagrangian (1B7|) . one can proceed 
to the Hamiltonian of the system (see Ref. [2] for the details of the derivation). 



H 



E 

i=l 



2fii 2 



dp 



^ 2 



+ 



c 



/ii + /g^ d[3v[3 



2r2[/xi(l - 0' + li2e + /o d(3u{(3 - C)1 /^i + /i2 + d(3u 

Extrema in the einbein fields are understood either in the Hamiltonian ( 1B8I) or, alternatively, 
in its spectrum. In the latter case the einbein field method is a variety of the celebrated 
variational method in Quantum Mechanics. 
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Now, if the contribution of the string to the total inertia of the rotating system (denomina- 
tor of the last, angular-momentum-dependent, term in the Hamiltonian flB8l) ) is neglected, 
then the extrema in all einbeins can be taken analytically yielding for the Hamiltonian (this 
procedure is exact for L = 0): 



H = y + ml + y fp + ml + ar, (B9) 
or, in the one-particle limit (mi = M ^ oo, m2 = m), 



where we omitted the infinite contribution of the static particle mass M. After a canonical 
quantisation of the Hamiltonian (IBlOp we reproduce the spinless Salpeter Eq. (^^. The 
role of the proper string dynamics in the Hamiltonian (]B8I) as well as supplying the spinless 
Hamiltonian with spin -dep endent terms are discussed in the literature — see, for example, 
Ref. 3l| and Refs. [s], [s^, respectively. Discussion of the proper string dynamics in the 
formalism of the Schwinger-Dyson-type Eq. ( fT2l) can be found in Ref. |17i]. Calculations 
of various hadronic spectra in the framework of the Hamiltonian (IBSP supplied by the per- 
turbative exchange and by the spin-dependent terms demonstrate a good accuracy of the 
predictions (see, for example, recent results for the spectrum of heavy-light D, Dg, B, and 
Bg mesons 3^). Notice that, in the case of light quarks, the major contribution to the 
spectrum of the Hamiltonian flB8l) comes from the confining QCD string. Therefore, this 
case can be referred to as the case of the "heavy" string (as opposed to the case of heavy 
quarks when the proper string dynamics gives only small corrections). 



